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S : Abstract 

^ I 74n explicit solution of the stationary one dimensional half-space boundary value 

O I problem for the linear Boltzmann equation is presented in the presence of an arbi- 

trarily high constant external field. The collision kernel is assumed to be separable, 
which is also known as "relaxation time approximation" ; the relaxation time may 
^ , depend on the electron velocity. Our method consists in a transformation of the 

0\ ', half-space problem into a nonnormal singular integral equation, which has an ex- 

r^ I plicit solution. 

o' 



I Introduction 

The motion of electrons or holes in a bulk semiconductor under the action of an 
external field E can be described by a probability density f{t,x,p) for finding a 
Q i charge carrier at time t and point x with momentum p. Under the assumptions of 

^ ' quasineutrality [Q] and a low electron concentration, so that interactions among the 

electrons themselves (via collisions or Pauli's principle) can be neglected 0, this 
^ ■ probability density satisfies the linear Boltzmann equation 



(dt+pd^- Edp)f{t,x,p) = Qfit,x,p) 



where we have set the charge and effective mass of the electrons equal to unity. The 
linear collision operator Q describes scattering of the electrons by vibrations of the 
host lattice and consequently strongly depends on the model for the electron-phonon 
interaction. 



In this article we consider collision models with separable collision kernels. These 
models are also known as relaxation time approximations with momentum depen- 
dent relaxation times. Although these models do not account for some high field 
effects in semiconductors , they are widely used because of their relatively simple 
mathematical structure. Whilst the solution of the initial value problem is easily 
obtained in these models Q 0, the boundary value problem is far more intricate. 
Up to now, a solution of the boundary value problem in the relaxation time ap- 
proximation has been found only by Cercignani for a specific p-dependence of the 
relaxation time ^ and by Dalitz in the zero temperature limit |0]. 

In this article, we present a solution of the stationary one dimensional half-space 
problem with a constant electric field for quite general p-dependent relaxation times. 
The somewhat simpler case of a constant relaxation time is treated in the unpub- 



lished preprint ||2J] in great detail. In the half-space problem the semiconductor is 
taken semiinfinite (x G (0, C)o)); at the boundary x = 0, electrons are shot in with a 
given momentum distribution (f. Thus we seek the solution of 

(1) {pd,-Edp)f{x,p) = Qf{x,p) 
with the boundary conditions 

(2) /(0,p) = v^(p) for p>0 and lim /(a;,p) = 

For zero electric field E = Case Q] has found a solution via an eigenfunction 
technique in the context of neutron diffusion in solids. Completeness theorems for 
these eigenfunctions have been proven by Zweifel [§ for the general equation 
h{p)dxf + Af = under the assumption that ^ be a selfadjoint positive operator 
in some appropriate function space. In the case of a nonzero electric field however, 
this general result is no longer valid because the operator A = —Edp — Q is not 
selfadjoint. 

Nevertheless for specific collision operators, an eigenfunction expansion of the 



solution might be possible also in the presence of a constant external field. In [10| a 



kind of half-range completeness of the eigenfunctions has been proven for a constant 



relaxation time in the zero temperature limit, and in [|T^ the same has been proven 
for two very specific collision kernels at arbitrary temperature of the semiconductor 
host lattice. In this article we show that the half-space problem in the relaxation 
time approximation generally can be solved via an eigenfunction expansion, provided 
the collision frequency is an analytic function of the momentum p that increases not 
more than linearly in p for p ^ oo. 

We shall prove that each boundary value (p that is a Laplace transform of any 
tempered distribution can be written as a superposition of "eigenfunctions" . This re- 
sult (see end of section V) is in complete analogy to Case's half-range completeness 
theorem for his "singular eigenfunctions" . Moreover we shall see that the eigenfunc- 
tion representation of the boundary value generally contains a singular contribution 
of the Maxwellian of the temperature of the semiconductor lattice. This is in agree- 



ment with a conjecture of Stichel and Strothmann |12 



II The collision model 

In the nondegenerate situation, the electron phonon colhsion operator Q in (1) has 
the general form 0] 

(3) Qfip) = M{p) J dp' K{p', p) f{p') - f{p) J dp' M{p') K{p, p') 

where M{p) := Mq e~^ /^^ is the Maxwellian at temperature of the semiconductor 
lattice, which is assumed to be in equilibrium. The collision kernel K{p,p') is a 
symmetric and positive distribution. The second integral in (3), which gives the 
total scattering rate, is called collision frequency z/; its reciprocal is called relaxation 
time T 

(4) uip) = :^--= jdp'M{p') K{p,p') 

Throughout this article we assume the collision kernels to be seperable, that is 
K{p,p') = u{j)) ■ u{p) = {t{p) ■ t{p'))^^. This assumption is also known as relaxation 
time approximation (RTA). If we normalize the Maxwellian so that / M{p) i>{p) dp = 
1, we obtain the RTA in its usual form 

(5) 2/(p) = ;^(M(p)/rfV^-/(p)) 

and by comparison with (4) we see that r(p) indeed is the relaxation time, because 
of our normalization J M/t = 1. The relaxation time approximation (5) generally 
can be a good approximation for the collision operator, provided the collision kernel 
K is a measurable function with a finite norm J dp J dp' M{p) M{p') K'^{p,p') < oo. 
For then (5) is equivalent to keeping only the largest term in the eigenfunction 
representation of the symmetric and square integrable kernel J M (p) K {p, p') J M {p') 
If 



III Eigenfunction expansion of the solution 

In order to solve the stationary Boltzmann equation (1), we make the separation 
ansatz fx{x,p) := e'^^^gx^p). We then obtain for the functions g\ the "eigenvalue" 
equation (Q + Edp) g\{p) = —^pg\{p), which reads with the collision operator (5) 

(6) [Edp + Xp- iy{p))gx{p) = -^{p) M{p) j dp' u{p') gx{p') 



In order to get rid of the integral on the right hand side of (6), let us normalize 
f V gx = 1- Then the general solution of (6) reads 

(7) gx{p) = ^ e-^^'/2^+^(^)/^ \cx - ] dqv{q) e(A~i^/e).^/2i?-^(.)/i^| 



where Cx is an integration constant and A^ (read "capital z/" ) denotes the primitive 
function of the colhsion frequency 

(8) N{p):= rdqv{q) 

Jo 

For general A, (7) is integrable with respect to p only if linip^oo ^{p)/p = c < oo. In 
fact, if this condition is violated, the only Li-integrable solutions of (6) are M{p) 
with A = E/6 and the homogeneous solution with A = 0. If the condition is satisfied, 
the solution (7) is Li-integrable for every A ■ sign(i?) > c. For simplicity we assume 
that c = 0. 

Because of the boundary condition at infinity (2), only positive values are allowed 
for A. Consequently, we must assume that E > 0, for otherwise g\ would grow 
exponentially for p -^ ±oo. This means that the field E must act in such a way 
that the electrons are driven back to the boundary x = 0. 

Inserting the eigenfunction (7) into the normalization condition Jdpu g\ = 1, we 
find for the constant Cx after a partial integration in the numerator 

7 (/pp e~^p'/^^+^(p)/^ / dq v{q) e^>^-Eie)q^liE-N{q)iE 

(9) Cx = 5^ '- 

/ dppe-^P^/'^^+^(py^ 

— oo 

Since equation (1) is linear, any superposition of solutions fx{x,p) is a solution of 
(1) too. Thus we assume the solution to be of the form 

/•oo /"OO 

(10) f{x,p)= d\A{\)fx{x,p)= d\A{\)e-'^gx{p) 

Jo Jo 

The expansion coefficients A{X) must be determined from the boundary value (/?. 
Setting a; = in (10) yields 

/•oo 

(11) f(0,p) = cp(p)= d\Ai\)gxip) for p>0 

Jo 

If this equation can be solved for A{X), then the half-space problem (1) and (2) is 
solved by (10). Hence the task is to determine the class of boundary values (/? for 
which (11) has a solution A{X) and to determine this solution in terms of ip. We shall 
see that (11) indeed can be solved if the boundary value is a Laplace transform of 
any tempered distribution, for then (11) is equivalent to a singular integral equation 
which can be solved explicitly. 

It is interesting to note that the representation (10) of the solution of the sta- 
tionary Boltzmann equation (1) can be written in a different form which is more 
general than the representation via eigenfunctions. If we define 

oo -. 

h{0 :=|rfAe-^«/^^A(A) (|"rfpz/(p)e-^^'/2^+^(^)/^y 



then f{x,p) = J^dXA{X) e gx{p) takes the form 



oo 



pN{p)/E 

f{x,p) = e^(P)/^ %2 ^ 2Ea;) + — ^ j dp' iy{p') e^'^P'^/'' x 



-oo 



Jdq u{q) e"'^(^)/^M(g) h{p^ + p'' - g' + 2Ex) 



The boundary condition at x = leads to an integral equation for the unknown 
function h. Maybe a solution of this equation is possible for a more general class of 
boundary values (p. 

IV A singular integral equation for the expansion 
coefficients 

If we insert the eigenfunctions (7) into the boundary condition f{0,p) = ^{p) for 
p>0 and make the substitutions t = p'^/2E and s = q^/2E, the integral equation (11) 
for the expansion coefficients A{\) reads 

E .„,^„„ ,_, 1 ,_, 



(12) A e-»<'^)/''^(v'2Bi) = fd\A{X)C^e- 



oo t 

- JdX A{X) e-^' jds e^'^E/2s u{V2E^) e-^(v^2^)-«^/^ 



Now let us make the assumption that the left hand side of (12) can be written as a 
Laplace transform of any tempered distribution \I^ 

(13) -^ e-^(^^)/^(^(v/2Et) = HdX e-^¥(A) 
Mq Jo 

Additionally let us assume that there is a continuous function Af with A/"(yu) > for 
/i>0 and 

(14) ^E/2s e'^^'^^'^/^iy{V2Ks) = /'"rf^e-^W(^) 

Examples for which this holds are z/(p) = ap" with — 1 < a < 1 and sums of these 
powers. If we integrate (14) from s to infinity, we find 

(15) E^-NiV2E-s)/E ^ dfie-^'W{fi)/fi 

Jo 

Therefore J\f/fi is locally integrable at /i = and consequently A/'(0) =0. Integrating 
a second time from s to infinity, we find lim^^oA/'//i = 0, so that /9(A) defined below 



(see 19) is a continuous function of A. Insertion of (13) and (14) into (12) and usage 
of the shifting theorem of the Laplace transform leads to (remember uniqueness of 
the Laplace transform) 



(16) ^(A) = {c> + p^. ,_^';'_^ }a(a) 

oo 

- e(A - E/e)j\f{\ - E/e) v fdfi 







/i — A 



where the symbol '"P" means Cauchy's principal value and denotes Heaviside's 
step function. This equation is a singular integral equation for the expansion coeffi- 
cients A{X). In the limit 6 —>■ 0, the G-function in front of the singular integral over 
A vanishes, which simplifies the situation considerably and allows the treatment 
described in 0, section 6. 

In contrast to "normal" singular equations however, the right hand side of (16) 
has a zero at A = E/6. Thus we are dealing with a "nonnormal" problem. This 
zero is easily split off by a partial integration using —EdqC"^/^ = ve~^/^ in the 
numerator of Cx (see 9) and usage of E = J^dfiAf/fi (see 15). The result is 



(17) ^{X) = (a - E/e) ■ I a{X) A{X) + /3(A) V fdfi 



/U — A 
with the notations 



;i8) a(A) = vfdfi 



^ /i(A - /i - E/e) 



\p^/2E+N{p)/E r J„ „ J\^E/e)q^/2E-Niq)/E 



J dppe-^P/'^+^^P>/'' Jdqqe 

-oo 



(19) p{X) = -Q{X-E/e) 



— oo 

Ar(A - E/e) 



x-E/e 
V Solution of the singular integral equation 

In appendix A.I we will prove that a(A) < for A < E/e. Therefore a{X) ± m[3{X) 
has no zeroes for A G IR"*", and we can transform the nonnormal problem (17) into 
a normal problem simply by dividing both sides by (A — Ee). 

However, if we divide the left hand side of (17) by (A — E/e), we obtain a 
distribution, even if ip' is a smooth function. Although the standard theory of singular 



integral equations |T^ only deals with Holder-continuous functions, physicists have 



for a long time been applying the same methods to distributions as well (see ||15 
for a comprehensive but heuristic treatment); meanwhile Estrada and Kanwal 
have presented a rigorous theory ^^. Recently the semiinfinite Hilbert transform 
of distributions has been discussed in a more general framework 



If we divide (17) by (A — E/6), we obtain the normal singular integral equation 



(20) T{X) = a{X)A{X) + l3{X)P dfi 



fi — X 



The tempered distribution T is the most general solution of the (distributional) 
equation T ■ (A — E/6) = \P. Thus we can write 

(21) T{X) = V^^^^+CE6{X-E/e) 

where ce is an arbitrary constant which will be essential later on. The symbol 
'"P" stands for "particular solution", which is known to exist from Hormander's 



theorem |]T^ . If we define principal value integrals in the more general sense of |T^ 
or [|l^], we may read this symbol as "principal value". 

The zero on the right hand side of (17) has introduced an arbitrary constant into 
our problem. This is in agreement with Prossdorf's result in a nondistributional 



framework ||T9| that such a zero increases the index of the equation by one. 



Although the problem of "nonnormality" is solved, there is still a problem with 
eq. (20). The standard theory of singular integral equations is only applicable to 
equations on a finite interval, but our equation lives on a semiinfinite interval. We 
can circumvent this difficulty with the transformation 

(22) /i = i^ and A = ^^ 

1 — y 1 — X 

which maps the interval 0< A<oo monotonously onto the interval — l<x< 1. This 

transformation does not seem to be well known. For instance Paveri-Fontana 

and ZwElFEL |^0[ recently presented an ab initio derivation of an inversion formula 



for the half-Hilbert transform which turns out to be equivalent to the well known 
inversion formula on a finite interval [^ if the transformation (22) is made. 
With this transformation, equation (20) transforms into 

(23) f(x) = a(x) A(x) + (3(x) V f dy ^^ 

y — X 



-1 



with the definitions 



Ti'-±^] A^' + ^ 



T{x):= ^^-^^ and A{x^ — ^^ ~ ^ 



X i — X 



a{x) := a ( — '■ — ) and P{x) := (3 



X J \l — X 



Although the theory of singular integral equations on a finite interval is well de- 
veloped (see [0 or p2| for a comprehensive summary), we will sketch the solution 
of (23) in some detail. Let us start with the definition of the analytic function F 
associated with the distribution A: 

1 }..>M^') 



F(z) := — f dx' ^^^ for zG(D\ [-1,11 
^ ^ 27riJ x' -z \ L > J 



which is a holomorphic function in (D\ [—1, 1]. If z approaches the cut [—1, 1] from 
above or below, its boundary values are given by the Plemelj formula 

1 } A(x') 1 

lim F(x±ie) = F^(x) = V dx' -^^ ± -A(x) 

^-0 ^ ^ ^ ^ 2m J x'-x 2 ^ ' 

Hence our integral equation (23) can be written as a relation between these two 
boundary values (remember that a ± Tii(3 ^ 0) 

(24) t^^^l— = mHimp^x) - F~{x) 

a{x) — 7ri(3{x) a{x) — 7ri(3{x) 

which is known as a Riemann-Hilbert problem in the literature. Now we are looking 
for a function x which is holomorphic in the complex plane cut at [—1, 1] and satisfies 
the boundary condition 

a{x) + 7ii(3{x) X^{^) 
a{x)-TTif3{x) X~{x) 

If we find such a function, we can determine x^ from (24) and therefore also F, and 
from A = F^ — F~ we can obtain A. A function with this property is 

(25) x{z) = (^ + 1)'"(2 - l)"e^(") with 

1 

(26) r{z) := - fdx'^^^ and 



71 J X — Z 

-1 



(27) 7(x) := —M\z ^ =arccot' ^' 



2i Va-vri/3/ Vvr/3(x)/ 

In (27) we are choosing the main branch of the arcus cotangens (see appendix A. II). 
The exponents m and n must be chosen in such a way that x ^i-nd F are integrable 



near the end points z = ±1 (see p2| for details; this is the crucial point where the 
interval must be finite): 

— 1 < m — liinx-*-i'~f{x)/7T <1 

— 1 < n + lima.^_|_i7(a;)/7r < 1 

8 



In appendix A. II it is proven that 7(— 1) = and 7(+l) = vr. Therefore we must 
choose m = and n = —1. Thus the index of the integral equation, which is defined 
hj m + n, is —1. Consequently, (23) has a solution A only if the inhomogeneity T 
satisfies the orthogonality relation [|l^] 

:28) [dx ^(^)^"(^) - 



a{x) — TTi(3{x) 

Fortunately T contains an arbitrary constant, which we will choose in such a way 
that this orthogonality relation holds. If (28) is satisfied, the solution of our integral 
equation (23) is unique and reads 

A{x) = F^{x)~F~{x) with 



Xiz) 27riJ^ a{x') - 7ri(3{x') {x' - z) 

By application of the Plemelj formula to (25)-(27) we can express x^ iii terms of 
principal value integrals. Moreover, the solution of the original equation on the 
semiinfinite interval (0, cx)) is easily obtained via the transformation inverse to (22). 
If we set X = {\ — 1)/(A + 1) and x' = (/i — l)/(/i + 1), we obtain 

/3(A) e-^W ^7^ (/x + 1) e^('^)r(/x) 



V jd^i 



with the notations 



(30) GiX) := -^Vfdfi 

n J 





(// + 1)(//-A) 



and 



31 g(a) := — In ^ = arccot — ^ 

The orthogonality relation for the inhomogeneity T becomes 

,32) U , ^'^'-°'" ^0 



According to Noether's theorem on singular integral equations p3|5 ^ necessary 



and sufficient condition for the solvability of the inhomogeneous equation (20) is the 



orthogonality of the inhomogeneity T to all solutions -B(A) of the adjoint homoge- 
neous equation 



oo 



(3{^^)B{^^) 



(33) a(A)i?(A)-p/rf^ ^^^^T =0 

J u — X 



u — A 
^ 

In our situation, this equation only has one solution. Comparison of Noether's 
orthogonality condition ^TB = o with (32) shows (the same result is obtained, of 
course, by direct solution of (33)) 

(34) B{X) = 

Because of , _a)'(a+i) ~ ~^ + a+T ^^'^ (32), we may write (29) as 

^(A).__i__L,)r(A)- mvU^Mm] 

Now let us insert T according to (21). The constant ce must be chosen 

1 7...o....^ n^ 



(35) ^---^7^/^/^^(^)^ 



so that the orthogonality relation (32) is satisfied. Finally we arrive at the solution 
of the singular integral equation for the expansion coefficients: 

(36) A{\) = ,,,, / ,,,,,, { a{X) V- '^^^^ 



a2(A) + 7r2/32(A) [ ' ' X-E/O 

where -B(A) is given by (34). 

Obviously, our solution f{x,p) = J d\A{X) e~^^g\{p) is well defined only if ^/a 
is locally integrable at A = 0. From eq. (A.l) in the appendix we conclude that for 
A^O 

l/a(A) ~ const ■ /" cippe-^^'/^E+^^P)/^ ^ oo 

J —oo 

Thus, if we demand 

oo oo 

(37) oo > /rfppv^(p)= /rfppv^(p)e-^(P)/^e^(P)/^ 





oo oo 



■\p^/2E+N{p)/E 



10 



the solution of the half-space problem is given by the eigenfunction expansion (10) 
with expansion coefficients (36). 

Let us repeat this result in the form of a half-range completeness theorem for 
the eigenf unctions g\. 

Theorem: The integral N over the collision frequency be such that 

Jo 

for some continuous function M with A/'(/i) > for fi > 0. Then each function 
if : IR^ —>■ ]R with f^ dp pip < cxo, that is a Laplace transform of any tempered 
distribution (p 

y,{p)= / dfi${fi)e-^P/^^ 
Jo 

is a superposition of eigenfunctions 

V{p)= / d\A{X)g^{p) 
Jo 

where the expansion coefficients are given by (36) with 

Remarks: a) Relation (38) between \I/ and <P is an immediate consequence of (14) 
and the convolution theorem of the Laplace transform. 

b) According to (36), the expansion coefficients A{X) generally contain a singular 
contribution from the eigenvalue A = E/9, the corresponding eigenfunction of which 
is the Maxwellian M{p) at temperature 9. This confirms a conjecture of Stichel 
and Strothmann |jl2[, who took this as a starting point for an asymptotic analysis 
of the boundary value problem. 

VI Examples 

VI. I Maxwellian as boundary value 

If the boundary value is a Maxwellian of the temperature 9 of the host medium, 
that is ip{p) = M{p) = Mqc"^^'^^, the solution of the half-space problem reads 
f{x,p) = e~^^^^M{p), as can be deduced from (1) with the use of QM = 0. Hence 
we should obtain A{X) = 6{X — E/9), so that this situation provides a nice test of 
our solution. 

Indeed, for (p = M the inverse Laplace transform W defined in (13) reads 

11 



Thus the expansion coefficients (36) read in this case 

a(A)/3(A)^ (3{X) ^ 7,„ 5(/i)/9(M) 



V fdfi 



X-E/e B{X) (A - E/9) y ^ fi-X 

and the constant ce is given by (note that the symbol V can be omitted because of 
/3(A) = for A < E/9) 



B{E/e)J •- p.-Eie 

Because of relation (33), the integrals over 5/3 can be replaced by a(A)-B(A) and 
a{Eie) B{E/e). If we take care oif3{E/e) = 0, we obtain indeed A{X) = 6{X-E/e). 

VI. II Constant relaxation time 

In case of a constant relaxation time z/(p) = const. = 1/r, the quantities a, (3 and 
B in the formula (36) for the expansion coefficients can be calculated explicitly. 
Moreover there is a relation between a and (3 that allows a simplification of some 
expressions. For details of the calculation we refer the interested reader to pl|. 



Without loss of generality we can set r and 6 equal unity |^ . Then the function 
M that is necessary for the calculation of ^ via (38) reads ^{ij) = {E/2TTfiY''^e~^''^^^ 
and the functions a and (3 are given by 



a{X) = J — X'/'e-'/''^'V fdx 

V TT J 



2X{X-E)-x 

[^ ~1/2E(X-E) 

where the principal value integral in a allows for an explicit evaluation 

~ a;-3/2g-i/x { yy^(l-./^e-y2J^e-''dt) for y < 



V dx ■ 

^ ~ ^ [ y\^ (l - '^^/ye-y life'' dt) for y > 

Moreover in this situation the function B{X), which generally contains a principal 
value integral (see (34) and (30)), can be expressed in terms of a regular integral. 
We find (apart from an irrelevant constant which cancels out in (36)) 

S(A) oc A-3/2eV2SA^(A) with 

gifi) 2X-E 



i{X) = exp 




fi + X~E 2fi-E 
12 



where g{fi) = aTccot{a/TTj3) is the function defined in (31). The function ^ is a 
smooth function of A which is very close to a straight fine; it is sketched in figure 1 
for different values of the field E. 

At first glance, the divergence of -B(A) for A -^ might cause trouble in the inte- 
grals in (36). However, it follows from (37) that \P{X)X~^'^e^'^^'^ is locally integrable 
at A = 0. Hence all integrals in (36) are well defined. 
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A Appendix 

A.I Negativity of a(A) for A < EjO 

For A < E jQ^ the first integral in the definition of a[X) (18) is a regular integral 
instead of a singular integral. Hence we may write this integral with the use of (15) 



oo 

— / dppe 



oo 

2 



{X~E/e)p^/2E~N{p)/E 

Thus, combining terms in (18), we obtain for A < E/6 

7 dpp e" V/2£+JV{p)/B Jdqq (>{>^~E/e)q^/iE-N{q)/E 
(A.l) a{X) = ---^ ^ ^ 

/ dppe-^p''/'^^+^(py^ 

— oo 

Since N{p) is an increasing function, the denominator of (A.l) is always positive. 
Now we will prove that the numerator is positive too. For this purpose, let us 
convert all integrals over regions with negative values of p and q into integrals over 
positive values. If we do so, we can write for the numerator of (A.l) 

oo oo 

{\-E/8)qy2E~N{q)/E 



Jdppe-'P"/^^+''^^/^Jdqqe 

P 

oo oo 

- Jdppe-'p"/'^+'"^-P^/^ Jdqqe 



p 

oo oo 

{X-E/0)q'^/2E-N(q)/E 





oo p 

Xp^/2E+Ni-^p)/E fj^ ^ JX~E/e)qy2E~N{~q)/E 



+ Jdppe-^P^'^+''^-P)^^ Jdqqe^ 



In the first term we write J°° dq = J^ dq — Jq dq. Collecting terms with /q°° dq and 
Jq dq, the numerator of (A.l) reads 



oo 



(A.2) jdppe-^^"'^'' Jdqqe'^^-^/'^""/'^ (e 



,{X-E/9)qy2E ( ^N{p)/E~N(q)/E _ ^N{~p)/E~N{q)/E 



> for p > 



p 

(X-E/9)qy2E f N(~p)/E-N(-q)/E _ M{p)/E-N{q)/E 



dqqe^ 



TTT' 

< tor p > g 



14 



We can estimate the positive first term in the curly braces via J^ dq > Jq dq. Thus 
we arrive at the following lower estimate for (A. 2): 

p 

^^-Xpi2E [^^^(\-EI0)q'/2E (N(-p)/E-N(-q)/E _ M(-p)/E^N{q)/E 



fdppe-^P'/^^ fdqqe^ 





> for g > 

Since this lower estimate is positive, the numerator in (A.l) is positive, and conse- 
quently a{\) is negative for A < E/9. 

A. II Behaviour of ex and (3 for large A 

The index of our singular integral equation depends crucially on the values of 

^ ' ^^ ' 2i \a{X) - 7ri(3{X) J 

for A -^ and A -^ oo. Thus we need to determine the behaviour of the argument 
of the logarithm as A varies from zero to infinity. 

First note that the absolute value of this argument is equal to unity for any 
value of A, because numerator and denominator are complex conjugates. Moreover, 
because oi f3 = for A < E/6, the argument is +1 for A < E/6 and with the choice 
ln(l) = we have 

^(A) = for A < E/e 

Now let us consider real and imaginary part of the argument of the logarithm 

, , ,. ^ fa + 7ii(3\ a^ — n'^P'^ , ^ ( a + -Kil3\ 27raf3 

(A.4) Re — = -^— — r-:- and Im 



a — nip j a^ + vr^/?^ \a — nijd j a^ + vr^/S^ 

We know from the previous section that a{E/6) is negative, and from the definition 
(19) we know that /3(A) < for A > E/6. In consequence the imaginary part is 
positive if X is slightly greater than E/6. 

Moreover, from (A.4) we see that the real part cannot be positive whilst the 
imaginary part is zero. In other words the positive real axis cannot be crossed for 
finite X. 

In order to determine the behaviour for A — ;► cxd, we use the lemma on the 
asymptotic behaviour of a and [3 (see below). An immediate consequence of the 
lemma is that limA->oo/3/a = and a{X) > for large values of A. Hence {a + 
Tiij3)/{a — 7rif3) approaches +1 with a negative imaginary part for A ^ oo. 

Applying this result to our function g{X) defined in (A. 3), we obtain the limiting 
values 

(A. 5) lim5f(A) = and lim g{X) = vr 
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The logarithm may be expressed in terms of the arcus tangens function 

f'^/3W\ ( a(A) \ 

(A. 6) q(X) = arctan — --— = arccot — -— — 

We must choose the branch of the arcus tangens, or arcus cotangens respectively, 

in such a way that (A. 6) is a continuous function of A and (A. 5) is satisfied. The 

resulting branches are sketched in figure 2. 

Lemma: As X approaches infinity, the asymptotic behaviour of a and (3 is given 

by 

(A.7) f3{X) ~ ■sJe/27tu{0)X-'^/^ 

— oo 

Proof: (A.7) is a consequence of (14) and the Tauberian theorems on the Laplace 
transform (see |2^, part 1), which may be applied because of A/" > 0. If we remember 
the definition f3{X + E/O) = Q{X)Af{X)/X, we see that (A.7) holds. 

The asymptotic evaluation of a{X) is a bit more intricate. Let us examine each 
term in the definition (18) separately: 

The first term on the right hand side of (18) behaves like 

(A.9) '- -^Jd^^^^^^ - ^ 



x-E/eJ "^ /x x~E/e 

The denominator of the second term in (18) can be written as an integral over the 
positive real axis, which is easily evaluated asymptotically with Laplace's method 
(see [^, chapter 3) 

oo 

(A. 10) /rfpe-"^'/2i=;p (^giv(p)/i? _ ^Ni-p)/E^^ a^oo ^(o)v/2^ A'^^/^ 
^ ' ' 

The numerator of the second term on the right hand side of (18) can be evaluated 
as follows: after conversion of all integrals over the negative real axis into integrals 
over the positive axis, we substitute t = q^/2E and s' = p'^/2E. If we then transform 
s = s' — t and change the order of integrations, we obtain 

oo oo 

E"^ fdse-^' jdte-'^^l^ [^N{yj2E{s+t))/E-N{V2Ei)/E _ ^N{-y/2Eis+t))/E-N{-V2Ei)/E 
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Again we may apply Laplace's method |]2^; we only need to determine the be- 



haviour of the inner integral for s ^ 0. The term in curly braces behaves like 



{...} -J' 



V 

s — 



V2Ei 



Thus the inner integral behaves like s J^^dpe^ ' i'{p)/E, and consequently the nu 
merator in (18) behaves like 



E 



oo 



(A.ll) '^^ ^ J dpe'^''/''u{p) 

— oo 

Collecting (A. 9) - (A.ll) and insertion into a{\) according to (18) yields (A. 8). 
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